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Abstract 


The  micro-Doppler  (m-D)  effect  is  caused  by  moving  parts  of  the  radar  target.  This 
effect  can  degrade  the  quality  of  the  ISAR/SAR  image.  However,  at  the  same  time, 
it  carries  information  about  the  features  of  moving  parts.  The  separation  of  the 
patterns  caused  by  the  stationary  parts  of  the  target  from  those  caused  by  the  mov¬ 
ing  (rotating  or  vibrating)  parts  is  the  topic  of  this  report.  Two  techniques  for  the 
separation  of  the  rigid  part  from  the  rotating  parts  have  been  proposed.  The  first 
technique  is  based  on  the  time-frequency  (TF)  representation  with  sliding  windows 
and  order  statistics  techniques.  The  first  step  in  this  technique  is  the  recognition  of 
rigid  parts  in  the  range/cross-range  plane.'  In  the  second  step,  the  spectrogram  and 
order  statistics  setup  are  employed  to  obtain  the  signals  caused  by  the  moving  parts. 
The  second  technique  based  on  the  Radon  transform  can  be  applied  in  the  case  of 
very  emphatic  m-D  effect.  In  the  first  step  the  rotating  parts  are  recognized,  based 
on  the  inverse  Radon  transform.  After  masking  these  patterns,  a  radar  image  with 
the  rigid  body  reflection  can  be  obtained.  The  proposed  methods  are  illustrated  by 
examples.  The  proposed  methods  not  only  focus  the  distorted  SAR/ISAR  images, 
but  also  provide  additional  information  about  the  rotating/vibrating  features  of  the 
target. 


Resume 


L’effet  micro-Doppler  (m-D)  est  cause  par  les  parties  mobiles  d’une  cible  radar.  II 
peut  degrader  l’image  ISAR/SAR.  Toutefois,  il  est  egalement  porteur  d’information 
sur  les  caracteristiques  de  ces  parties  mobiles.  La  separation  des  formes  dues  res- 
pectivement  aux  parties  fixes  et  aux  parties  mobiles  (toumantes  ou  vibrantes)  de  la 
cible  est  I’objet  du  present  rapport.  Deux  techniques  ont  ete  proposees  pour  separer 
ces  deux  types  de  formes.  L’une  utilise  la  representation  temps- frequence  (TF)  avec 
fenetres  mobiles  et  statistiques  d’ordre.  Elle  permet  d’abord  de  reconnaitre  les  par¬ 
ties  rigides  dans  le  plan  site/azimut.  Dans  un  deuxieme  temps,  le  spectrogramme  et 
la  configuration  des  statistiques  d’ordre  sont  employes  pour  obtenir  les  signaux  dus 
aux  parties  mobiles.  La  deuxieme  technique,  basee  sur  la  transformee  de  Radon, 
peut  etre  appliquee  dans  le  cas  d’effets  m-D  tres  marques.  Les  parties  toumantes 
sont  d’abord  reconnues,  a  l’aide  de  la  transformee  de  Radon  inverse.  Apres  mas- 
quage  de  ces  formes,  on  peut  obtenir  une  image  radar  du  corps  rigide.  Des  exemples 
illustrent  les  methodes  proposees,  qui  non  seulement  permettent  de  mettre  au  point 
les  images  SAR/ISAR  deformees,  mais  foumissent  une  information  supplemental 
sur  les  parties  toumantes/vibrantes  de  la  cible. 
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Executive  summary 


If  a  target  or  any  structure  on  a  target  has  mechanical  vibration  or  rotation  in  addi¬ 
tion  to  its  bulk  translation,  it  might  induce  a  frequency  modulation  on  the  returned 
signal  that  generates  sidebands  about  the  target’s  Doppler  frequency  shift.  This  is 
called  the  micro-Doppler  effect.  Radar  signals  returned  from  a  target  that  incor¬ 
porates  vibrating  or  rotating  structures,  such  as  propellers  of  a  fixed-wing  aircraft, 
rotors  of  a  helicopter,  or  the  engine  compressor  and  blade  assemblies  of  a  jet  air¬ 
craft,  contain  micro-Doppler  characteristics  related  to  these  structures.  The  micro- 
Doppler  effect  enables  us  to  determine  the  dynamic  properties  of  the  target  and  it 
offers  a  new  approach  for  the  analysis  of  target  signatures.  Micro-Doppler  features 
serve  as  additional  features  that  are  complementary  to  those  made  available  by  ex¬ 
isting  methods.  The  micro-Doppler  effect  can  be  used  to  identify  specific  types  of 
vehicle,  and  determine  their  movement  and  the  speed  of  their  engines.  Vibrations 
generated  by  a  vehicle  engine  can  be  detected  by  radar  signals  returned  from  the 
surface  of  the  vehicle.  From  micro-Doppler  modulations  in  the  engine  vibration 
signal,  one  can  distinguish  whether  it  is  a  gas  turbine  engine  of  a  tank  or  the  diesel 
engine  of  a  bus. 

Micro-Doppler  effect  appears  in  the  SAR/ISAR  image  of  a  target  whenever  the 
target  has  one  or  more  rotating  or  vibrating  parts.  If  the  frequency  modulations 
on  the  returned  signal  caused  by  the  moving  parts  are  not  filtered,  then  the  micro- 
Doppler  effect  can  introduce  distortion  in  the  SAR/ISAR  images.  The  observation 
of  very  large  distortions  from  experimental  SAR/ISAR  data  has  been  reported.  On 
the  other  hand,  the  micro-Doppler  effect  also  carries  information  about  the  features 
of  moving  parts  of  a  target  that  are  complementary  to  existing  target  recognition 
methods.  This  report  proposes  two  techniques  for  the  separation  of  rigid  body 
radar  image  from  micro-Doppler  effect  caused  signal  in  SAR/ISAR  imaging.  The 
first  technique  is  based  on  the  order  statistics  of  the  spectrogram  samples.  In  order 
to  separate  the  patterns  caused  by  the  rotating  and  rigid  body  parts,  the  spectrogram 
is  evaluated  for  various  window  sizes.  The  spectrogram  evaluated  with  a  narrow 
window  produces  a  high  concentration  of  the  signal  caused  by  rotating  parts,  while 
the  spectrogram  evaluated  for  a  wide  window  produces  a  high  concentration  of  the 
signal  corresponding  to  the  rigid  body  parts.  The  second  technique  is  based  on 
the  Radon  transform  processing  of  the  radar  signals.  It  can  be  useful  in  the  case 
of  radar  targets  with  a  dominant  micro-Doppler  signals.  Both  proposed  methods 
have  shown  satisfactory  accuracy  for  the  considered  complex  simulated  examples, 
including  rotating  and  vibrating  target  parts. 

The  proposed  method  not  only  can  focus  distorted  SAR/ISAR  images,  but  also  can 
provide  additional  information  about  the  rotating/vibrating  features  of  the  target. 
Therefore,  new  algorithms  and  methods  will  have  to  be  further  investigated  to  sep- 
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arate  the  stationary  and  non-stationary  components  of  the  target  of  interest. 


T.  Thayaparan;  2006;  Separation  of  target  rigid  body  and 
micro-Doppler  effects  in  ISAR/SAR  imaging;  DRDC  Ottawa 
TM  2006-187;  Defence  R&D  Canada  -  Ottawa. 


IV 


DRDC  Ottawa  TM  2006-187 


Sommaire 


Lorsqu’une  cible  ou  une  structure  d’une  cible  presente  une  vibration  ou  rotation 
mecanique  en  plus  d’une  translation  globale,  elle  peut  induire  une  modulation  de  la 
frequence  du  signal  reflechi,  ce  qui  produit  des  bandes  laterales  de  part  et  d’autre 
du  deplacement  de  frequence  Doppler  de  la  cible.  On  parle  alors  d’effet  micro- 
Doppler.  Les  signaux  radar  reflechis  par  une  cible  comportant  des  structures  vi- 
brantes  ou  toumantes,  par  exemple  helices  d’aeronef  a  voilure  fixe,  rotors  d’helico- 
pteres  ou  ensembles  compresseur-aubes  d’avions  a  reaction,  presentent  des  carac¬ 
teristiques  micro-Doppler  associees  a  ces  structures.  L’effet  micro-Doppler  permet 
de  determiner  les  proprietes  dynamiques  de  la  cible  et  offre  une  nouvelle  approche 
pour  l’analyse  de  signatures  de  cibles.  Ces  caracteristiques  completent  celles  que 
foumissent  les  methodes  existantes.  L’effet  micro-Doppler  peut  servir  a  identifier 
des  types  particuliers  de  vehicules  et  a  determiner  leur  mouvement  et  la  vitesse  de 
leurs  moteurs.  Les  vibrations  produites  par  le  moteur  d’un  vehicule  peuvent  etre 
detectees  par  des  signaux  radars  reflechis  a  la  surface  du  vehicule.  A  partir  de  mo¬ 
dulations  micro-Doppler  du  signal  de  vibration  du  moteur,  on  peut  distinguer  s’il 
s’agit  du  moteur  a  turbine  a  gaz  d’un  tank  ou  du  moteur  diesel  d’un  autobus. 

L’effet  micro-Doppler  apparait  dans  l’image  SAR/ISAR  d’une  cible  lorsque  celle- 
ci  comporte  une  ou  plusieurs  parties  toumantes  ou  vibrantes.  Si  les  modulations 
de  frequence  du  signal  reflechi  dues  aux  parties  mobiles  ne  sont  pas  filtrees,  l’ef¬ 
fet  micro  Doppler  peut  deformer  les  images  SAR/ISAR.  On  a  observe  des  defor¬ 
mations  tres  importantes  de  donnees  SAR/ISAR  experimentales.  Par  contre,  I’effet 
micro-Doppler  foumit  aussi  de  l’information  sur  les  caracteristiques  des  parties  mo¬ 
biles  d’une  cible  qui  complete  l’information  provenant  des  methodes  de  reconnais¬ 
sance  de  cibles  existantes.  Le  present  rapport  propose  deux  techniques  pour  separer 
1’image  radar  du  corps  rigide  et  le  signal  dfl  a  l’effet  micro-Doppler  dans  l’image- 
rie  SAR/ISAR.  La  premiere  se  fonde  sur  les  statistiques  d’ordre  d’echantillons  de 
spectrogramme.  Pour  separer  les  formes  dues  aux  parties  toumantes  et  aux  parties 
rigides,  le  spectrogramme  est  evalue  avec  diverses  tailles  de  fenetres,  une  fenetre 
etroite  produisant  une  forte  concentration  du  signal  due  aux  parties  toumantes,  et 
une  fenetre  large  une  forte  concentration  du  signal  correspondant  aux  parties  du 
corps  rigide.  La  deuxieme  technique  se  fonde  sur  le  traitement  des  signaux  radar  par 
la  transformee  de  Radon.  Elle  peut  etre  utile  dans  le  cas  de  cibles  radar  avec  signal 
micro-Doppler  dominant.  Les  deux  methodes  proposees  ont  montre  une  exactitude 
satisfaisante  dans  les  cas  de  simulation  complexe  etudies,  comprenant  des  parties 
de  cible  toumantes  et  vibrantes. 

La  methode  proposee  peut  non  seulement  mettre  au  point  les  images  SAR/ISAR 
deformees,  mais  aussi  foumir  une  information  supplemental  sur  les  parties  tour- 
nantes  et  vibrantes  de  la  cible.  De  nouveaux  algorithmes  et  methodes  devront  done 
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etre  etudies  plus  en  profondeur  afin  de  separer  les  composantes  fixes  et  mobiles  de 
4a  cible  d’interet. 
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1  Introduction 


Micro-Doppler  (m-D)  effects  appear  in  the  SAR/ISAR  image  when  the  target  has 
one  or  more  rotating  or  vibrating  parts.  This  effect  can  degrade  the  quality  of  the 
radar  image.  However,  the  m-D  effect,  at  the  same  time,  carries  information  about 
the  features  of  moving  parts  (type,  velocity,  size,  etc.).  Several  papers  have  been 
written  about  the  ways  to  deal  with  the  m-D  effect.  The  wavelet  analysis  of  the  heli¬ 
copter  and  human  data,  along  with  the  time-frequency  (TF)  representation  based 
imaging  system,  is  presented  in  [1,  2,  3].  Details  on  the  physics  of  the  m-D  effect, 
with  some  typical  examples,  are  given  in  [4,  5, 6,  7].  A  method  for  separating  the  m- 
D  effect  from  the  radar  image,  based  on  the  chirplet  transform,  is  proposed  in  [8,  6]. 
The  part  of  the  radar  signal,  obtained  as  a  result  of  the  m-D  effect,  after  appropriate 
separation  from  the  rest  of  the  radar  target  image,  can  be  used  for  estimating  the 
rotating  part  parameters.  In  [8,  6],  the  received  signal  has  been  expanded  into  a  set 
of  chirplet  functions.  The  obtained  chirplet  functions  are  divided  into  two  groups: 
the  chirplets  with  the  small  chirp-rate  parameter  (associated  to  the  rigid  part)  and 
the  chirplets  with  the  large  chirp-rate  parameter  (associated  to  the  moving  parts). 
The  calculation  burden  of  this  procedure  is  high,  since  the  chirplet  dictionary  could 
be  extremely  large.  Both  wavelet-based  and  chirplet  based  procedures  are  used 
in  [9,  6]  to  extract  the  m-D  features  such  as  the  rotating  frequency  of  the  antenna 
from  SAR  data.  The  analysis  of  the  TF  representation  application  in  radar  target 
identification  is  presented  in  [10].  The  reduced  interference  distributions  from  the 
Cohen  class  are  applied  as  a  tool  for  the  target  identification.  A  technique  for  the 
m-D  effect  estimation  from  the  reflected  signal,  based  on  the  TF  signatures  and  de¬ 
composition  of  basis  functions,  is  presented  in  [1 1],  This  technique  can  be  used  for 
m-D  effect  signals  that  can  be  represented  as  sinusoidal  FM  signals. 

In  this  report,  two  procedures  for  the  separation  of  signals  coming  from  the  rigid 
(slowly  moving)  body  and  rotating  parts  are  proposed.  In  order  to  analyze  the  m-D 
effect,  we  will  perform  TF  (more  accurately  space/spatial-frequency)  analysis  using 
a  sliding  window.  Based  on  the  order  statistics  of  the  obtained  TF  representations, 
we  will  make  decisions  on  whether  the  components  belong  to  the  rigid  body  or  to 
the  rotating  target  parts  [21],  The  other  approach  is  based  on  the  Radon  transform 
(RT).  The  TF  representation  of  the  sinusoidal  FM  signal  can  be  concentrated  in  a 
single  point  by  applying  the  inverse  RT  to  the  TF  representation.  In  this  way,  the 
patterns  of  the  rotating  parts  can  be  easily  extracted  from  the  target  signal  [21]. 

The  report  is  organized  as  follows.  ISAR  signal  modeling  is  described  in  Section  2, 
modeling  signals  from  both  the  rigid  body  and  the  moving  parts  of  the  target.  The 
method  for  separation  of  the  rigid  body  and  the  moving  elements  from  the  radar 
image,  based  on  the  order  statistics  of  the  TF  representation  with  a  sliding  window, 
is  given  in  Section  3.  The  separation  of  rotating  patterns  from  the  radar  image  using 
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the  inverse  RT  is  presented  in  Section  4.  The  simulation  study  is  given  in  Section 
5.  The  concluding  remarks  with  possible  directions  for  the  further  research  are 
given  in  Section  6.  Annex  A  provides  the  basic  procedure  of  the  order  statistics 
estimations. 
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2  Radar  Signal  Model 


Consider  a  radar  signal  consisting  of  M  continuous  wave  coherent  pulses: 


M—l 

vM(t)  =  ^2v0(t-mTr),  (1) 

771=0 

where  v0(t)  is  the  basic  impulse  limited  within  the  interval  —Tr/2  <  t  <  Tr/2. 
The  linear  frequency  modulated  (FM)  signal  is  used  in  our  simulations  as  a  basic 
impulse:  vq (t)  =  exp(j-KBt2 /Tr),  where  B  is  the  bandwidth  control  parameter 
while  Tt  is  the  pulse  repetition  time.  The  alternative  radar  model  used  in  practice 
has  radar  pulses  with  stepped  frequencies,  however,  this  model  is  not  studied  in  this 
report.  The  defocusing  effect  considered  in  this  report  and  time-frequency  (TF) 
signatures  of  the  obtained  radar  signals  have  similar  behavior  for  these  two  forms 
of  radar  signals  [  1 9,  20], 

The  signal  emitted  toward  the  radar  target  can  be  written  as: 


u(f)  =  e?2vIolvM{t),  (2) 

where  /0  is  the  radar  operating  frequency.  The  received  signal,  reflected  from  a 
single  reflector  target  at  distance  d(t),  is  delayed  for  2d{t)/c,  with  cr  being  the 
reflection  coefficient  and  c  being  the  propagation  rate: 


Ufi(t )  =  au(t  —  2 d(t)/c) .  (3) 

The  demodulation  of  the  received  signal  can  be  performed  by  multiplying  the  re¬ 
ceived  with  the  transmitted  signal  u(t): 

q(t )  =  au*(t  —  2  d(t)/c)u(t)  = 


M-l  M- 1 

a exp(j'47r fod(t) / c)  ^  Vg(t  -  2 d{t)/c  -  mTr )  ^  v0(t  -  mTr  -  T0).  (4) 

m  =  0  777=0 

The  parameter  T0  is  used  in  radar  imaging  for  compensation  of  target  distance.  For 
properly  selected  T0  and  after  highpass  filtering,  the  signal  q(t)  can  be  approxim¬ 
ately  written  as: 


DRDC  Ottawa  TM  2006-187 


3 


A-/-1 

q(t)  ~  <rexp(j4irfod(t)/c)  ^  v${t  —  2d(<)/c 

m= 0 


A/-1 


m'Tr)vo(t  -  mTr)  =  E  q{m,t), 

(5) 


jn.=0 


where 


q(m,  t)  =  a  exp(^47r fod(t) J c)vq (t  —  2 d(t)/c  —  mTr)vo(t  —  mTr), 
t  e  [(m  —  1/2 )Tr,  (m  4-  1/2 )Tr), 

=  crexp(j47r/o<i(i)/c)  exp(j4irBd(t)(t  —  mTr)/[cTr))  exp(  —  iirB{2d(t)/c)2 /Tr). 

(6) 

Keeping  in  mind  that  B  -C  /o,  we  can  neglect  exp(—  jir B (2d(t) / c)2 /Tr)  with  re¬ 
spect  to  the  other  two  components.  The  value  of  q(m,t)  can  approximately  be 
written  as: 


q{m,t)  «  <jexp(j4'!rfod(t)/c)exp(j4iTBd(t)(t  —  mTr)/(cTr)).  (7) 

This  signal  is  commonly  given  in  the  form: 


q(m,T)  «  cr  exp(j4ir fod(r  +  mTr)/c)  exp(j4nBd(r  +  mTr)T / (cTr)) ,  (8) 

where  t  =  r  +  mTr.  Parameter  r  €  [— Tr/2,Tr/2)  is  referred  to  as  fast-time, 
while  m  =  0, 1, M  —  1,  is  called  the  slow-time  coordinate.  Commonly,  in  actual 
radar  systems,  signals  are  discretized  in  the  fast-time  coordinate  with  sampling  rate 
T,  =  Tr/N,  t  =  nTs,  where  n  E  [— 1V/2,  N/ 2).  The  classical  radar  setup  assumes 
that  the  radar  target  position  is  a  linear  function  of  time  d(t)  —  D0  +  Vt.  Then  the 
radar  model  produces: 


q(m, t )  ~  (7  exp(j47r/o A)/ c)  exp(j4TrVrnf0Tr/c )  exp{j4-rrTBD0/{cTr)).  (9) 


A  two-dimensional  (2D)  FT  of  this  signal  over  m  and  r  is  approximately: 


Q{u> 


r  M-\ 

T,Um)  =  /  ^2 

Jr  rn= 0 


re 


-jbjTT-ju)mrn 


dr 
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-A/- 1 

/  ^  <7  exp(j'47r/oiDo/c]  exp{j4nVmfoTr/c)  exp(j47rr B Do/ (cTr))e~JUTT~:,WmmdT 

JT  m= o 

.  A/-1 

=  cr  exp(j47r/oDo/c)  /  exp(j4-KTBDo/(cTT))e~JUJrrdT  exp(j'47r\/m/oTr/c) 

'^T  771=0 

N  /  •  ,  f  n  /  u/  ,  nn  //  ~^Vh%/c)M/2) 

=  (27r)(7 exp(j/47r/o J90/ c)8(ut  -  ^BDq/{cTt))— — j -  .  .  w  . 

sin((u;m  *±kv  jol?/ 0)1*2) 


?~JWni  111 


Xe~~j(um-4nVf0Tr/c)(M-l)/2 


(10) 


\Q(u)T,Um)\  =  \{2TT)a  exp(j4nf0D0/c)S(ujT  -  ^-kBDq/ (cTr)) 


x 


sin((^m  47rK/pTr/c)M/2)  j^m_4nV  fnTr/c)(M-i)/2 
sm((u>m  -  AttV foTr/c)/ 2) 


=  |27r||cr||  exp(j4TrfoDo/c)\\d(ur-4irBD0/(cTr))\ 


sin((u;m  -  ^ttV f0Tr/c)M/2) 
sin((wm  -  4nV  f0Tr/c)/2) 


x\e-j(u,m-4nVhTr/cW-m 


(11) 


=  |  (27T  )g5(ut  —  47T  BDq/(cTt)) 


sin((t<;m  -  47rVr/0Tr/c)M/2) 
sin((wm  -  4ttV f0Tr/c)/2) 
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«  (27r)cr<S(wr  -  AirBD0/(cTr))M8(ujm  —  AnV /0Tr/c ) 


since 


sin(aM) 

sin(a) 


M8{a) 


for  relatively  large  M. 

Therefore,  for  large  M  we  can  write  the  magnitude  of  Q(u>T,u>m)  as: 
|Q(wt,wto)|  ~  (2ir)a8(uiT  -  AirBD0/{cTr))M8(ujm  -  2 Vf0Tr/c). 


(12) 


The  distance  can  be  approximately  written  as  d(t)  ~  R(t)+xp  cos (9(t))+yp  sin 
where  R(t)  is  the  distance  of  the  target  rotation  center  from  the  radar,  where  co¬ 
ordinates  of  the  scatterer,  for  r  =  0,  are  (x p,yp)  (see  Figure  1).  The  coordinate 
system  is  formed  in  such  a  way  that  the  coordinate  x  is  the  line  of  sight.  Assume 
constant  rotation  velocity  9(t)  =  uRt,  with  relatively  small  angular  movement  of 
the  target  \uRTr\  -C  1  (it  implies  that  cos (9(t))  ~  1  and  sin {9(t))  as  0).  Accord¬ 
ing  to  the  introduced  conditions  d(t)  ~  xp  and  v(t)  =  d'(t)  =  —xp9'(t)  sin {9(t.))  + 
yp9'(t)  cos (0(f))  w  yp9'(t)  cos {9(t))  «  ypuiR.  Commonly,  it  is  assumed  that  R(t)  is 
compensated  by  adjusting  T0  in  (4).  Thus,  we  will  not  consider  it  in  our  algorithm. 
Then  \Q{uT,ujm)\  can  be  written  as: 


|Q(u;r,u;m)|  «  (2ir)aM5(ojT  -  Air Bxp/ (cTr))8{ujm  -  Airypu Rf0Tr / c) 

=  ( 2n)aM8{ojT  -  c\xp)8(ujm  -  C2yp).  (13) 

It  represents  the  ISAR  image  of  scatterer  ( xv ,  yp)  for  a  given  instance  under  in¬ 
troduced  assumptions.  Note  that  the  constants  that  determine  the  resolution  of  the 
radar  image  are  given  by  C\  =  AnB/(cTr )  and  C2  =  ATru}RfoTr/c.  The  radar  image 
is  formed  as  superposition  of  radar  images  of  all  scatterers  (xp,  yp),  p  =  1,  2, ...,  P. 
It  is  approximately  given  as: 


p 

| i^(cUr ,  CUm) |  —  ^](27r)crp5(wr  C\Xp)8((AJm  —  c2yp),  ()^) 

P=1 


where  crp  is  the  reflection  coefficient  that  corresponds  to  the  p-th  scatterer  point. 

The  common  ISAR  imaging  model  assumes  that  all  the  point-scatterers  share  the 
same  angular  motion.  However,  the  radar  target  can  have  moving  parts  (rotating  or 
vibrating).  In  this  model,  the  received  signal,  within  the  m-th  radar  sweep,  can  be 
written  as  [8]: 
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r  s  ^  j-  x 

uo  =  E^  exp  l  —j — —  [i?B(f)  +  xpcos 9B(t)  +  yps\ndB(t)]  >  + 
p=i  L  c  J 

o r  exp  l  -j—^-[RR{t)  +  xRcosdR(t)  +  yRsin9R(t)]  >  , 

where  ap  is  the  reflection  coefficient  of  the  p-th  scatterer,  /  is  the  radar  operat¬ 
ing  frequency,  (xp,  yp)  are  coordinates  of  the  scatterer  point  (£  coordinate  is  in  the 
direction  of  line-of-sight,  according  to  Figure  1).  The  target  translation  and  angu¬ 
lar  motion  are  denoted  by  R(t)  and  0(t),  with  indices  B  and  R  corresponding  to 
the  rigid  body  and  rotating  parts,  respectively.  Here  we  assume  a  single  rotating 
point.  In  general,  we  will  have  a  radar  target  with  an  arbitrary  number  of  moving 
parts.  The  motion  compensation  techniques  [8,  12]  are  employed  to  remove  the 
influence  of  the  translation  motion.  Thus,  we  can  assume  that  RB(t)  — ►  0  and 
RR(t)  — >  0.  For  the  main  body  holds  \9B(t)\  <C  1.  It  results  in  cos 9B(t)  ~  1 
and  sin  9B(t)  ~  9B(t )  =  ujBt,  where  ujb  is  the  effective  body  rotation  rate,  after 
the  motion  compensation.  This  approximation  cannot  be  applied  to  rotating  parts 
since,  dR(t )  can  vary  rapidly.  The  received  signal  can  now  be  written  as: 


u\  f  -^Tr/o 

Qm{t)  =  2^^Pexp  <  -J  — 

P=i  ^  c 


Xp  +  ypuBt] 


4tt/0  . 


-{-o' R  exp  -j — —  [xrcos  dR(t)  +  yRsindR(t)] 
c 


(16) 


The  first  term  in  (16)  represents  sinusoids  concentrated  at  frequencies  proportional 
to  yp,  p  €  [1,  P].  The  information  on  other  coordinates  can  be  obtained  based  on 
processing  of  all  returned  radar  signals.  The  radar  image  is  commonly  obtained 
by  a  2D  FT  of  the  compressed  received  signal  denoted  as  q(m,t )  =  qm(t).  The 
second  term  in  (16)  represents  sinusoidally  modulated  FM  signal  caused  by  the 
rotating  part  with  a  relatively  large  0R(t)  during  the  repetition  interval.  This  signal 
can  cover  the  large  portion  of  the  TF  plane.  It  can  mask  a  significant  part  of  the 
radar  image.  The  influence  of  the  rotating  and  vibrating  parts  on  the  radar  image  is 
commonly  referred  to  as  the  m-D  effect.  The  detailed  considerations  of  physics  of 
this  effect,  mathematical  derivations  and  simulation  studies  are  given  in  [4,  5,  7], 
This  effect  should  be  removed  from  the  radar  image  of  the  rigid  body  part,  but 
at  the  same  time  it  contains  important  information  related  to  the  radar  target.  In 
the  next  section,  a  technique  for  the  separation  of  the  rigid  body  and  rotating  parts 
influence  from  the  received  signal  is  presented.  It  is  based  on  the  running  window 
TF  representation  and  the  L-statistics.  In  Section  4  the  RT  is  applied  for  removing 
influence  of  strong  rotating  parts. 
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3  Time-frequency  analysis  and  order 

statistics  approach _ 

3.1  Basic  idea 

In  order  to  explain  the  approach  for  separation  of  the  m-D  effect  from  the  ISAR 
image,  we  will  consider  the  experimental  setup  used  in  [8].  The  target  with  five 
scatterers  at  positions:  (0,0),  (±6,0)  and  (0,  ±6)  is  considered  (positions  are  in 
meters).  The  reflection  coefficient  for  the  central  point  is  5,  for  the  second  pair  of 
points  is  1,  while  for  the  third  pair  is  2.  The  positions  of  the  six  scatterers  are  shown 
in  Figure  2.  Scatterer  6  rotates  around  scatterer  2  at  a  rate  of  6.67  Hz  and  a  rotation 
radius  of  20  cm,  with  the  reflection  coefficient  of  3.33.  We  assume  that  the  radar 
has  a  centre  frequency  /  =  10  GHz,  bandwidth  B  =  800  MHz,  and  PRF  =  1400 
Hz.  The  target  body  rotates  about  4°  over  384  pulses  during  the  data  collection  time. 
Therefore  the  coherent  integration  time  is  Tr  =  0.2743  s.  The  range  and  cross-range 
resolutions  are  Rrange  =  c)2B  =  0.1874  m  and  RCross-range  =  c/2  fucTc  =  0.2149 
m,  respectively. 

For  the  analysis  of  the  received  signals  q(r,  rn)  we  will  use  the  spectrogram  (SPEC) 
as  a  TF  tool: 


SPEC(t.ojt',  rn) 


q{r  +  Tw,m)iv{Tw)e 


dr.,, 


(17) 


where  w(tw)  is  a  window  function.  We  will  assume  that  the  signal  q(T,m)  is 
discretized,  producing  discretized  SPEC{n ,  cvr;  m). 

The  SPEC  for  various  window  lengths  exhibits  different  properties  for  various  kinds 
of  signals.  The  SPEC  of  the  considered  radar  signal,  for  a  narrow  Hanning  window 
with  32  samples,  is  shown  in  Figure  3a.  The  SPEC  calculated  with  a  wide  Hanning 
window  with  128  samples  is  depicted  in  Figure  3b.  In  both  figures  we  can  see  the 
stationary  patterns,  caused  by  the  rigid  body,  parallel  to  the  time-axis.  The  SPEC 
of  the  rigid  body  parts  is  better  concentrated  for  wider  than  for  narrower  windows. 
Also,  a  sinusoidal  pattern  caused  by  the  rotating  point  can  easily  be  seen.  It  is  in 
accordance  with  the  assumed  model  for  these  parts.  However,  the  sinusoidal  FM 
signal  is  better  concentrated  for  the  narrower  window  case. 

This  example  clearly  shows  that  the  rigid  body  and  rotating  parts  in  the  TF  plane 
should  be  studied  with  the  TF  representation  with  various  window  lengths.  Upper 
arrow  in  Figure  3a  denotes  the  frequency  where  only  the  influence  of  the  rotating 
parts  exists.  For  that  frequency  relatively  small  number  of  TF  samples  assume  high 
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Figure  2:  Point-scatterer  representation  of  the  original  target  consisting  of  five  rigid 
points  (1-5)  and  one  rotating  point  (6)  with  strengths  2,  5,  2,  1,  1,  and  3.33, 
respectively. 


values.  However,  in  the  region  of  the  signal  caused  by  the  rigid  body  (see  lower 
arrow)  all  TF  samples  have  high  values. 

For  a  narrow  window  in  the  SPEC  we  can  see  that  the  number  of  points  belonging 
to  the  rotating  parts  is  not  high  as  compared  to  the  considered  interval.  Therefore, 
we  can  consider  the  mean  value  of  the  smallest  SPEC  samples,  as  a  good  measure 
of  the  rigid  body  and  rotating  parts  influence.  In  the  region  of  the  rigid  body  this 
measure  will  be  high,  since  the  rigid  body  has  influence  over  the  entire  considered 
interval  while  the  influence  of  the  signal  caused  by  rotating  part  will  be  removed. 


3.2  L-statistics 

For  this  kind  of  analysis  a  linear  combination  of  order  statistics  (L-statistics)  gives 
a  powerful  notational  and  analytic  tool.  The  L-statistics  of  the  SPEC  samples  can 
be  defined  as  [13]: 


/V-l 

^  ^  rri)  (18) 

n=0 

where  S(n)(u/r;rn)  are  samples  from  the  set  SWrim  =  {SPEC(t,ut\  m),  r  6 
[0,Tr)},  sorted  into  a  non-decreasing  sequence:  m)  <  S(„+i)(wT;  rn).  In 

order  to  visualize  how  this  measure  can  be  used  to  distinguish  signals  caused  by 
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Figure  3:  Spectrogram  of  the  radar  signal  reflected  from  the  target  with  emphatic 
m-D  effect  calculated  with:  (a)  Narrow  window;  (b)  Wide  window. 
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Figure  4:  TF  representation  and  order  statistics  of  the  received  signal:  (a) 
Spectrogram;  (b)  L\{uit  );  (c)  L‘2 (u^ t  ) ’ (d)  L&  )  ■  The  x-axis  corresponds  to  time 
and  the  y-axis  corresponds  to  frequency. 


various  parts  of  the  radar  target,  we  repeat  Figure  4a.  Three  order  statistics  of  this 
SPEC  are  presented  on  the  right  hand  side  of  this  figure  (Figures  4b,  c  and  d).  They 
are  calculated  as: 

•  an  —  4/N  for  n  €  [0,  N/4)  and  an  =  0  elsewhere,  i.e.,  for  the  smallest  mag¬ 
nitude  SPEC  samples  (statistics  L\  k)); 

•  a/v/2-i  =  aN/2  =  1/2  and  an  =  0  elsewhere,  i.e.,  medians  of  the  SPEC  values 

for  fixed  uT  (statistics  and 

•  an  =  4/N  for  n  €  [3iV/4,  (V)  and  an  =  0  elsewhere,  i.e.,  for  the  highest 
magnitude  SPEC  samples  (statistics  L3(u;T)). 

It  can  be  seen  that  for  the  median  (L2(u;T)),  Figure  4c,  the  frequency  region  of  sta¬ 
tionary  components  is  noticeable.  However,  from  this  illustration  it  is  difficult  to 
recognize  that  there  are  three  components  along  the  y-direction.  From  Figure  4b, 
again,  the  region  of  stationary  components  can  be  seen.  In  addition,  from  three 
peeks  we  can  make  the  conclusion  that  scatterers  are  located  at  three  different  y- 
positions.  In  addition  to  the  complex  sinusoidal  components  in  the  third  case,  Fig¬ 
ure  4d,  the  sinusoidal  FM  signal  region  can  be  seen.  These  three  simple  functions 
will  be  used  in  the  construction  of  indicators  for  the  region  of  rigid  body  and  the 
region  of  rotating  parts.  Several  modifications  will  also  be  made  in  order  to  take 
the  advantage  of  the  two-dimensional  structure  of  the  signal  g(r,  m),  as  well  as  to 
keep  the  calculation  burden  within  a  reasonable  limit. 
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3.3  Algorithm 

Let  us  introduce  the  following  notation:  I{uiTruJm)  =  Iriq(ujt, 
where  I(ojt,  ujm)  is  a  function  indicating  the  region  in  the  range/cross-range  domain 
with  components  caused  by  target  (both  from  the  rigid  and  rotating  parts),  while 
4[G(wnWm)  is  the  region  of  rigid  part;  and  /Mov(wnwm)  is  the  region  of  moving 
parts.  In  general,  /rigO^t,  ^m)  0  /MOv(wr,u;m)  ^  0  holds.  The  domain  (ujT,ujm) 
will  be  referred  hereafter  as  the  range/cross-range  domain,  since  it  can  be  easily 
transferred  to  the  positions  (x,  y)  of  the  radar  target. 

As  an  indicator  of  target  in  the  range/cross-range  domain  the  following  simple  cri¬ 
terion  will  be  used: 


I(u>T,Um)  =  {  \  >  A  IN  IQ(Wr,Wm)|  (19) 

v  '  \  0  elsewhere, 

where  Q(uiT,u>m)  is  a  2D  FT  of  q(r,  m).  This  is  a  very  conservative  measure,  since 
the  global  threshold  used  here  is  equal  to  the  mean  of  \Q(ur,ujm)\,  calculated  in 
the  entire  range/cross-range  domain. 

Consider  a  2D  SPEC  of  the  received  signal: 


SPEC(r,uiT;  m,ujm;  Nw,  Mw)  = 


A/-1 

H  9(r  + 

Am=0 


tw,  m  +  Am)w(rw.  Am)e-jWrTw-jUmArndTw 


2 


(20) 


The  parameters  Nw,  Mw  are  added  to  denote  the  size  of  the  window  function 
w(Tw,Am)  in  the  corresponding  directions.  The  window  width  Mw  is  selected 
as  Mw  ta  M  (number  of  chirps).  However,  the  narrow  windows  in  the  r  direction 
are  used  in  order  to  control  the  concentration  of  the  sinusoidal  FM  patterns  caused 
by  the  rotating  points.  Namely,  if  the  window  is  very  wide  the  TF  representation 
of  the  sinusoidal  FM  signal  will  look  like  several  quasi  stationary  components  [14], 
see  Figure  3b.  By  decreasing  the  window  width  we  are  achieving  a  better  concentra¬ 
tion  of  these  components  around  the  instantaneous  frequency  of  the  sinusoidal  FM 
components.  Of  course,  if  we  decrease  the  window  length  too  much,  the  frequency 
resolution  will  be  decreased  and  semiperiods  of  this  signal  will  be  connected,  res¬ 
ulting  again  in  TF  representation  degradation. 

The  above  analysis  confirms  that,  in  determination  of  the  indicator  function,  the 
SPECs  should  be  calculated  with  various  window  lengths.  We  will  assume  that 
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window  lengths  are  selected  from  a  set  Nw  €  N,  for  constant  Mw.  The  SPEC  will 
not  be  evaluated  for  each  (r,  m),  but  only  for  equidistantly  spaced  instants  along  the 
r  direction.  The  positions  are  separated  for  A Tnw  (for  different  windows  different 
intervals  are  used;  in  our  calculation  study  they  are  A T^w  =  TCN/NW,  Nw  <  N). 
Then,  for  each  point  in  the  range/cross-range  domain  we  calculate: 


N- 1 

Li,Nw  (wr,  U>m)  ^  ^  ^n^(n)  ?  ^m)>  2  (21) 

i=0 

where  S(„)(uT,ujm)  are  elements  from  the  set: 

SNw  -  { SPEC(r,ujT;m,um ;  NW,MW),  r  =  r'  +  IATNw,1  €  Z }  (22) 

sorted  into  a  non-decreasing  order  5(n)(u;r,u;m)  <  5(7[+1)(u/T,u;m).  The  index  i 
denotes  an  earlier  introduced  form  of  L-statistics,  while  Nw  denotes  length  of  the 
window  function.  Based  on  the  above  considerations,  we  get  the  function  for  the 
rigid  part  (sinusoidal  components): 


'RIG 


(^rj^rn)  1 1  ^ Niu  ^m)i 


Nw  CN 


(23) 


where  /'Vm  (coT ,  cvm )  can  be  defined  as 


(^r  j  m  ) 


(wT,tUm)  >  I  AND 

1  £fc=-fC|  l^i,fVu)(t«'T+fcAu;T,wm+tAu;rn)|  (2K,i  +  l)(2L1 +1) 

E^-k,  |ii,w«(«r+fcAwT^m+/Awra)|  “  (2Ar2  +  l)(2L2+l) 

0  otherwise. 

(24) 


The  first  condition  in  (24)  is  equivalent  to  (19),  while  the  second  condition  is  in¬ 
troduced  to  further  emphasize  the  local  maxima,  corresponding  to  the  scattering 
points.  It  compares  values  of  Li,/v,„  samples  (smallest  SPEC  samples)  in  a  close 
neighborhood  of  the  considered  point  (uT,  ujm)  with  those  in  a  wider  neighborhood 
(tfi  <  K2  and  Li  <  l2). 

For  the  widest  window  from  the  set,  the  function  \Li^w(uT+k\  Au>T,ojm+l\  Aium,)| 
would  be  well  concentrated  around  the  frequency  of  the  rigid  body  signal  compon¬ 
ents.  The  sinusoidal  FM  components,  that  are  spread  in  the  TF  plane  (see  Fig¬ 
ure  3a),  could  have  high  values  of  this  function.  However,  for  a  narrow  window, 
the  function  \L\ ,nw(ut  4-  k\Au>T,ojm  +  /jAo;m)|  would  be  spread  around  the  rigid 
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body  components,  while  the  rotation  components  will  be  eliminated.  Therefore, 
the  product  of  I'Nw{u)T, cuin)  in  (23)  will  produce  the  function  that  is  equal  to  1  for 
a  very  narrow  region  in  the  range/cross-range  domain  (as  in  the  case  of  wide  win¬ 
dow)  with  eliminated  influence  of  the  rotation  parts  (as  in  the  case  of  using  narrow 
window). 

The  remaining  problem  is  determining  the  patterns  representing  the  moving  objects. 
The  indication  of  moving  object  parts  is  described  as: 


(25) 


where: 


/" 


N„ 


(cvr ,  com) 


1  ^3,/Viy  i  ^  (1  T  P)  {^t  i  ^m) 

0  otherwise. 


(26) 


The  indicator  function  fj(u(a;r,a;m)  will  be  equal  to  1  for  regions  in  the  range/cross¬ 
range  domain  where  the  ratio  between  the  highest  and  the  median  SPEC  samples  is 
higher  than  (1  +  /?),  where  p  is  a  threshold  value.  In  our  analysis  P  =  1  is  used. 
This  indicates  that,  along  with  the  stationary  pattern,  there  could  exist  a  pattern 
caused  by  the  moving  objects,  increasing  the  values  of  the  SPEC  at  some  instants 
of  the  TF  plane.  The  remaining  problem  is  to  remove  possible  influence  of  the  side 
lobes  caused  by  the  stationary  pattern  from  the  product  /(wT,wm)/f(iii(wT,wm).  To 
this  aim  we  form  a  function  i\(ojT,ujrn)  as: 


A-(cU-7-,  OJjrp  —  1  dllate(/R|G (uir,iOm))  (27) 

where  “dilate”  denotes  the  dilatation,  a  morphological  operation  that  extends  ar¬ 
gument  matrix  (under  extension  we  assume  that  the  neighboring  region  around 
^RiG(wx,wm)  =  1  becomes  also  equal  to  1)  [15].  Function  dilate(/R[C(wr,cvm)) 
for  3x3  closest  samples  in  the  range/cross-range  domain  can  be  defined  as: 

i  i 

di late (/RiG (wT,a;m))  \J  /rig(cut  T  -F  (28) 

k=-\ <=-i 

where  V/L-i  a(^)  =  a(~  1)  V  a(0)  V  a(l).  This  form  is  a  dilatation  with  struc¬ 
tural  element  square.  We  have  introduced  the  dilatation  since  the  indicator  func¬ 
tion  /RiG(wr,u;m)  is  very  sharp.  It  also  includes  corresponding  peaks  only  in  the 
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Step  1. 

Calculation  of  SPEC(T,uT;m,com ;  Nw,  Mw),  r  =  r'  +  lATNmtl  €  Z  in 
equidistantly  spaced  points  for  windows  Nw  from  set  N  that  contains  both 
wide  and  narrow  windows.  The  interval  A7/v,„  is  selected  to  be  A Tnw  = 
TrN/Nw  where  N  is  the  widest  window  from  the  set. 

Step  2. 

Determination  of  the  radar  image  as  2D  DFT  of  signal  q(r,  rn),  Q(cjt,  ojm). 

Step  3. 

Determination  of  region-of-nterest  J(wT,wm)  by  (19). 

Step  4. 

Calculation  of  order  statistics  Li  ^w(uiT,ujm),  i  =  1,2,3,  (21),  and  for  all 
the  windows  from  set  N. 

Determination  of  the  indicator  functions  I'Nm(ujT,ojm)  by  using  (24). 

Step  6. 

Calculation  of  region-of-interest  function  for  rigid  body  /RiG(wT,u;m)  = 
/(cur,u;m)  ri/v,„cN 

Calculation  of  (ujT,uim)  by  using  (26). 

b-'TTOM 

Calculation  of  A (ujt,  ojm)  by  using  (27). 

Step  9. 

Evaluation  of  the  indicator  function  for  the  part  of  range/cross¬ 
range  domain  caused  by  moving  parts  /Mov(wT,wm)  = 

/(wr ,  U>m) I (wT,  Wm)  A(wr ,  dm) ■ 

Step  10. 

Signal  caused  by  rigid  body  and  moving  parts  can  be  reconstructed 

as:  qRiG{t,m)  =  IFT{Q{u>T,u!m)IMG{ujT,(jjm)}  and  qMov(t,m)  = 

I FT ,  UJjTfi) I MOV  (^r ,  CJjn)  }  . 

Table  1:  Procedure  for  the  radar  image  analysis. 


range/cross-range-domain.  Therefore,  the  function  A (u>T,u>m)  allows  that  points 
that  are  not  close  to  the  rigid  body  part  could  be  identified  as  rotating  parts. 

For  the  separation  of  the  rotating  object  component  we  can  use  the  product 

Also,  as  it  can  be  seen,  in  the  case  of  moving  parts  de¬ 
tection,  we  did  not  employ  the  multiwindow  approach,  but  only  a  single,  relatively 
narrow,  window.  Namely,  in  this  case,  results  obtained  for  extremely  narrow  and 
wide  windows  are  not  useful  in  the  separation  process,  while  for  other  windows 
results  are  very  similar. 

Note  that  the  algorithm  is  derived  based  on  the  assumption  that  the  rigid  body  influ¬ 
enced  signal  can  be  modeled  as  a  sum  of  complex  sinusoids.  For  more  complicated 
target  motion  patterns,  the  motion  compensation  techniques  for  focusing  the  distor¬ 
ted  data  should  be  applied  in  the  first  stage  of  the  algorithm  [16]. 

The  procedure  for  separation  is  summarized  in  Table  1. 
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4  The  Radon  transform  based  algorithm  for 
the  separation  of  sinusoidal  FM  signals 


The  second  approach  that  will  be  presented  in  this  section  is  based  on  the  RT.  It 
can  be  used  only  for  the  rotating  parts  of  the  target  that  produce  pure  sinusoidal  FM 
signals.  This  approach  can  be  applied  only  in  the  case  when  the  influence  of  the 
rotating  parts  is  emphatic. 

The  RT  of  a  point  with  position  (xr,  yr)  is 


Rr(p,0)  =  CrS{p  -  Arsin{9  -  ifT)),  (29) 

where  cT  corresponds  to  the  reflection  of  the  point,  A  —  \J +  yf  and  <pr  = 
tan ~i{yr/xT).  The  period  of  the  RT  in  9  is  2?r.  It  means  that  the  inverse  RT  of 
a  sinusoidal  pattern  in  the  (9,p)  plane  is  a  point  in  the  ( xT,yT )  domain.  Since  the 
rotating  parts  produce  sinusoidal  patterns  in  the  TF  domain,  it  means  that  the  inverse 
RT  of  the  TF  representation  of  the  received  signal  with  m-D  effects  will  have  peaks 
at  the  positions  corresponding  to  the  rotating  parts.  Note  that  for  the  inverse  RT  we 
should  know  the  period  of  the  sinusoidal  pattern,  which  is  not  known  in  advance. 
In  addition,  the  radar  signal  may  consist  of  several  rotating  parts  with  different 
unknown  rotation  frequencies.  To  solve  this  problem  we  assigned  the  considered 
time  interval  (in  the  TF  representation)  to  a  set  of  possible  periods  in  9.  The  whole 
TF  representation  is  considered  as  periodic  in  time  with  Tp  corresponding  to  periods 
in  9,  for  example  from  7t/3  up  to  4tt  with  a  step  of  O.Itt.  The  value  of  47r  means 
that  exactly  two  periods  of  sinusoidal  pattern  exist  in  the  TF  representation. 

When  the  assumed  period  matches  the  true  one,  the  peak  in  the  inverse  RT  is  de¬ 
tected  at  a  position  defined  by  A  =  \JrAt  +  y%  and  <pr  =  tan ~l{yr/xr).  In  this 
way,  all  the  parameters  of  the  rotating  part  are  known,  including  the  frequency  of 
the  rotation  that  corresponds  to  the  period  assigned  for  the  case  where  the  peak  is 
found.  If  there  are  more  than  one  rotating  part,  with  different  or  the  same  rotation 
speed,  more  peaks  will  be  detected.  By  filtering  the  region  around  the  peak  and 
calculating  the  RT  of  the  filtered  inverse  RT,  we  get  a  clear  TF  representation  of  the 
signal  resulting  from  the  rotating  parts. 
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5  Results 


5.1  Example  1: 

In  the  first  example  we  considered  the  simulation  setup  described  in  Section  3  and 
used  in  [8],  The  radar  image,  including  the  notable  m-D  effect,  is  depicted  in  Figure 
5a.  The  algorithm  described  in  Section  3  is  applied.  The  Hanning  windows  are 
employed  along  both  directions  with  Mw  =  M  and  the  set  of  window  lengths  in 
N  =  {512,  256,  128,  64,  32,  16,  8}  is  used.  In  the  calculation  of  /rig(^t>  wm),  the 
expression  (23)  is  applied,  where  the  product  FI nw c n  ^nw  (wt ,  )  is  calculated 
iteratively  starting  from  the  widest  toward  narrower  windows  from  the  set  N.  If  the 
resulting  product  is  not  changed  in  consecutive  iterations  we  can  stop  the  procedure. 
This  is  an  important  step,  since  it  significantly  reduces  the  calculations  in  the  case  of 
targets  without  m-D  effect.  The  obtained  radar  image  of  the  rigid  body  is  depicted  in 
Figure  5b.  The  parameters  K\  =  L\  =  1  and  K2  =  L2  =  3  are  used  in  (24).  It  can 
be  seen  that  excellent  accuracy  is  achieved.  We  calculated  the  region  of  sinusoidal 
FM  component  by  using  (25),  with  Nw  =  64.  The  region  of  sinusoidal  FM  signals 
is  given  in  Figure  5c.  In  order  to  demonstrate  the  accuracy  of  this  approach,  the 
reconstruction  of  the  sinusoidal  FM  signal  for  m  =  0  is  given  in  Figure  5d.  The 
sinusoidal  FM  component  can  easily  be  recognized  from  this  illustration.  At  the 
same  time  some  important  features  of  the  moving  object  can  be  extracted. 

5.1.1  Vibration 

This  setup  has  been  used  for  the  alternative  experiment.  The  same  target,  but 
without  rotating  parts,  has  been  considered.  However,  the  entire  target  oscillates 
in  the  x  direction  with  amplitude  of  0.1  cm  and  frequency  of  oscillations  30  Hz. 
The  obtained  results  are  presented  in  Figure  6.  A  significant  improvement  in  radar 
imagining  is  achieved  with  the  proposed  technique.  At  the  same  time,  this  tech¬ 
nique  can  be  used  for  extraction  of  the  vibration  parameters.  Note  that  an  analysis 
of  m-D  effect  for  vibrating  targets  in  SAR  systems  is  given  in  [9],  [17]. 


5.2  Example  2 

A  more  complex  object,  with  45  scatterers:  {(k,  0)  for  k  e  [—6,6],  (0 ,k)  for  k  € 
[—6,  —1]  U  [1,6],  (k,±(6-k))  for  k  €  [1,5]  and  (k,  ±(6  +  k))  for  k  e  [-5,-1]}, 

is  considered.  All  scatterers  have  the  same  reflection  coefficient,  set  to  2.  There  are 
two  rotating  parts.  The  first  rotating  part  is  the  same  as  in  the  previous  example, 
while  the  second  rotating  part  rotates  at  a  rate  of  2  Hz  and  a  rotation  radius  of  0.5  m 
from  center.  The  reflection  of  the  second  rotating  part  has  been  set  to  7.  The  same 
analysis,  with  the  same  algorithm  setup,  like  in  the  previous  example  is  performed. 
Obtained  results  are  depicted  in  Figure  7.  It  can  be  seen  that  the  reconstruction  of 
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Figure  5:  Separation  of  rigid  body  and  m-D  effect  for  simple  case:  (a)  Q(u>r.  um); 
(b)  Reconstructed  radar  image;  (c)  Region  of  support  for  the  sinusoidal  FM  signal; 
(d)  SPEC  of  the  reconstructed  m-D  signal. 


objects  has  been  done  in  a  very  accurate  way.  At  the  same  time,  the  reconstruction 
of  the  signal  caused  by  the  rotating  parts  is  good. 


5.3  Radon  transform  based  algorithm 

In  this  case  we  also  apply  the  inverse  RT,  assuming  that  the  TF  representation  pat¬ 
tern  is  periodic  with  periods  from  7r/3  to  37t,  with  step  of  tt/18.  Two  peaks  of 
the  inverse  RT  concentration  are  detected  at  two  different  periods,  amplitudes  and 
phases,  see  Figure  8.  Calculating  the  RT  of  the  inverse  RT  filtered  around  these 
peaks,  we  obtained  the  representation  presented  in  Figure  9,  producing  very  accur¬ 
ate  reconstruction  of  the  rotating  parts  in  the  TF  domain. 
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Figure  6:  Separation  of  rigid  body  and  m-D  effect  caused  by  object  vibrations:  (a) 
Q(ojT,um);  (b)  Reconstructed  radar  image;  (c)  Region  of  support  for  the  signal 
caused  by  vibrations;  (d)  Spectrogram  of  the  reconstructed  m-D  signal. 
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Figure  7:  Separation  of  rigid  body  and  m-D  effect  for  two  moving  points:  (a) 
Q(ur,ujm);  (b)  Reconstructed  radar  image;  (c)  Region  of  support  for  the  sinusoidal 
FM  signal;  (d)  SPEC  of  the  reconstructed  m-D  signal. 
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30 


Figure  8:  Concentration  measure  of  the  inverse  RT  of  the  radar  signal  TF 
representation  as  a  function  of  the  assumed  angle  number  (top),  with  the  inverse 
RT  at  two  positions  where  maxima  are  found  (bottom). 


Figure  9:  TF  representation  corresponding  to  the  rotating  parts  in  the  last  example, 
obtained  by  using  filtered  inverse  RT.  The  x-axis  corresponds  to  range  time  and  the 
y-axis  corresponds  to  frequency. 
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5.4  UH-1D  Helicopter 


In  this  example  we  consider  the  simulated  signal  of  a  German  Air  Force  Bell  UH- 
1D  Helicopter  known  also  as  ‘Iroquois’.  The  simulation  is  performed  according 
to  [4],  Several  effects  are  emphasized  in  the  TF  representation  Figure  10.  The 
stationary  patterns  along  the  time-axis  correspond  to  the  rigid  body  reflection,  the 
vibration  of  the  target  or  the  radar-clutter  caused  by  the  moving  in  the  target  back¬ 
ground.  The  motion  of  the  two  main  blades  is  modeled  by  two  rotating  reflectors, 
producing  sinusoidal  FM  signals  with  a  large  magnitude  in  the  frequency  direction, 
(15).  The  main  rotor  flashes  are  simulated  by  the  signals  producing  lines  connect¬ 
ing  peaks  of  the  sinusoidal  FM  signal  with  time  axis.  The  smaller  pulses  that  can  be 
seen  in  the  right  hand  side  of  Figure  10a  correspond  to  the  tail  rotor  flashes.  These 
flashes  correspond  to  periodic  alignment  of  the  main  and  tail  rotors  to  maximally 
reflect  the  radar  signal.  Note  that  other  effects  that  can  be  observed  in  a  radar  image, 
including  multipath,  are  not  considered  here. 

Therefore,  the  simplified  model  of  the  reflected  UH-1D  signal  can  be  written  as: 


x(t)  =  Xfuc(t)  +  Xftorit)  +  xpi_ni(t)  +  xpi_r{t), 

where  x R/c(f),  XRor(t),  xpi_M(t)  and  Xpi_-r(t)  represent  signals  caused  by  the 
rigid  body,  rotation  of  the  main  rotor,  and  the  main  and  tail  rotor  flashes,  respect¬ 
ively.  The  signal  is  considered  within  the  interval  of  400  ms,  sampled  with  a  rate 
of  A t  =  1/48000  s.  Four  sinusoidal  components,  caused  by  the  rigid  body,  are  at 
the  frequencies  —10.3  kHz,  —2.5  kHz,  2.3  kHz  and  2.7  kHz.  Two  components  at 
—0.4  kHz  and  0.4  kHz  correspond  to  the  modulated  time  tones  commonly  added  to 
the  data  tape  [18].  The  sinusoidal  FM  signals,  corresponding  to  the  rotation  of  the 
main  rotor  blades,  are  modeled  as: 


XRor{t)  =  <7ROT[exp(j27rAfiorsin(27rf/Tflor))+exp(-j27rAfiorsin(27rf/Tfior))]- 

(30) 

where  Trot  =  175  ms  and  Arot  =  19  kHz.  The  main  and  tail  rotor  flashes  are 
modeled  as  broadband  pulses  given  as: 


XFL_M{t )  =  opi  m  ]^<5(t  —  kTRor/ 2)  *T  ^fl_a/(£) 

k 

XFL_r{t)  =  a FL_T^2,fi(t  ~  kTrAiL/2)  *r  hpi  rit),  (31) 

fc 
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where  T^ail  in  our  experiment  is  Ttail  =  35.8  ms,  while  liFL_M(t )  and  hpurit) 
are  cut-off  filters  given  in  the  frequency  domain  as: 


1  |cj|  <  2 TtAfiOT 
0  elsewhere, 


Hfl_t{u) 


1  2tt(7.35  kHz)  <  u  <  2tt(15.7  kHz) 
0  elsewhere. 


The  signal  is  corrupted  with  the  moderate  Gaussian  noise.  To  compare  our  simu¬ 
lated  radar  image  with  the  real  one,  refer  to  [4,  18]. 

The  previously  used  form  of  the  L-statistics  can  be  applied  for  separation  of  the 
rigid  body  and  rotating  parts  influence.  In  the  case  when  we  want  to  separate  the 
different  effects  caused  by  the  rotating  parts,  the  modified  L-statistics  can  be  ap¬ 
plied.  The  L-statistics  corresponding  to  the  smallest  SPEC  samples  is  presented  in 
Figure  10b.  As  a  threshold  we  adopted  10%  of  the  maximal  value  of  the  corres¬ 
ponding  L-statistics.  The  dotted  line  represents  the  region  of  detected  stationary 
patterns.  For  the  detection  of  the  tail  blades  region  we  have  considered  L-statistics 
of  30%  the  highest  SPEC  samples  for  each  frequency,  excluding  10%  the  highest 
ones  of  them  (Figure  10c).  Again,  the  region  is  estimated  based  on  the  values  that 
are  above  1 0%  of  the  maximal  value  of  the  L-statistics,  but  with  excluded  stationary 
components  regions.  Finally,  the  highest  10%  of  SPEC  samples  are  given  in  Figure 
lOd.  Based  on  this  illustration,  it  is  easy  to  determine  the  region  of  the  main  blades 
(comparing  to  the  threshold  determined  as  in  previous  cases,  but  not  excluding  any 
sample).  The  separation  of  the  target  components  is  shown  in  Figure  1 1. 

Note  that  the  proposed  approach  can  be  used  for  helicopter  radar  signals  where 
relatively  long  FM  signal  components  appear.  In  the  case  when,  in  the  recording, 
short  burst  components  dominate,  on  high  frequencies,  without  a  sinusoidal  FM 
component,  wavelet  transform  analysis  setup  can  be  efficiently  used  [1], 
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Figure  10:  Radar  imaging  of  the  simulated  helicopter  signal:  (a)  TF  representation; 
(b)  Detection  of  the  stationary  patterns;  (c)  Detection  of  the  tail-blades;  (d) 
Detection  of  the  main  blades  region. 
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Figure  11:  Sorted  spectrogram  samples  to  obtain  L-statistics 
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Figure  12:  Separation  of  components  of  the  target:  stationary  signal  pattern  (top, 
left);  main  rotor  flashes  (top,  right);  tail  rotor  flashes  (bottom,  left);  rotating  blades 
(bottom,  right). 
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6  CONCLUSION 


Two  techniques  for  separating  the  rigid  body  radar  image  from  micro-Doppler  ef¬ 
fect  caused  signal  in  ISAR/SAR  imaging  are  proposed.  The  first  technique  is  based 
on  the  order  statistics  of  the  SPEC  samples.  In  order  to  separate  the  patterns  caused 
by  the  rotating  and  rigid  body  parts,  a  spectrogram  is  evaluated  for  various  window 
sizes.  The  spectrogram  evaluated  with  a  narrow  window  produces  a  high  concen¬ 
tration  of  the  signal  caused  by  rotating  parts,  while  a  spectrogram  evaluated  for  a 
wide  window  produces  a  high  concentration  of  the  signal  corresponding  to  the  rigid 
body  parts.  The  second  approach  is  based  on  the  Radon  transform  processing  of  ob¬ 
tained  radar  signals.  It  can  be  useful  in  the  case  of  the  radar  targets  with  a  dominant 
micro-Doppler  effect  caused  signals.  Both  proposed  methods  have  shown  satisfact¬ 
ory  accuracy  for  the  considered  complex  simulated  examples,  including  rotating 
and  vibrating  target  parts. 

This  report  provides  the  preliminary  ground  work  for  this  challenging  field  of  re¬ 
search.  There  is  room  for  improvement  of  the  algorithm  in  terms  of  accuracy  and 
adaptability.  Two  components  of  the  algorithm,  i.e.,  spectrogram  and  sorting  of 
samples,  can  be  further  studied  in  order  to  improve  the  performance  of  the  al¬ 
gorithm.  For  example,  fast  sorting  procedures  such  as  quicksort  or  insertion  sort 
can  be  investigated  and  employed  to  reduce  the  calculation  burden  of  the  proced¬ 
ure. 

This  report  also  demonstrates  that  whenever  the  target  has  one  or  more  rotating 
or  vibrating  parts,  the  frequency  modulations  on  the  returned  signal  can  introduce 
large  distortion  in  the  SAR/ISAR  images.  On  the  other  hand,  the  m-D  effect  also 
carries  information  about  the  features  of  moving  parts  of  the  target  that  are  com¬ 
plementary  to  those  existing  recognition  methods.  The  proposed  method  not  only 
can  focus  the  distorted  SAR/ISAR  images,  but  also  can  provide  additional  inform¬ 
ation  about  the  rotating/vibrating  features  of  the  target.  Therefore,  new  algorithms 
and  methods  will  have  to  be  further  investigated  to  separate  the  stationary  and  non¬ 
stationary  components  of  the  target  of  interest. 
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Annex  A 
L-statistic 


This  section  provides  a  step  by  step  procedure  of  the  L-statistics  method  described 
in  Section  3  to  separate  the  rigid  and  rotating  parts  of  the  target.  The  procedure 
is  given  here  in  numerical  example  to  visualize  the  methodology.  Assume  that 
we  have  a  signal  with  three  sinusoidal  components  and  a  sinusoidal  FM  signal. 
Here,  we  will  consider  the  synthetic  time-frequency  plane  with  just  21x21  samples. 
We  will  also  assume  that  the  spectrograms  calculated  in  the  tables  are  with  wide, 
middle  and  narrow  size  windows.  The  goal  is  to  separate  the  sinusoidal  components 
(caused  by  the  rigid  body)  and  sinusoidal  FM  signal  (caused  by  a  rotating  part). 

Consider  a  simplified  matrix  of  the  signal  values  in  the  time-frequency  given  in  the 
next  table.  As  it  will  be  shown  later  the  same  situation  appears  in  practical  real 
examples. 
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The  blue  color  represents  sinusoidal  component  caused  by  rigid  body.  The  red 
color  is  for  the  part  caused  by  sinusoidal  FM  signal,  here  distorted  (due  to  the  wide 
window)  to  quasi-sinusoids.  The  green  color  are  the  samples  influenced  by  both 
sinusoidal  FM  signals  and  sinusoids.  Assume  that  for  each  frequency  we  perform 
the  sorting  in  time  from  the  smallest  toward  the  highest  values.  In  this  case,  since 
the  values  in  time  do  not  change,  we  obtained  in  fact  the  same  time-frequency 
representation. 
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Sorted  time  frequency  representation 
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The  smallest  samples  are  used  to  indicate  the  region  of  the  sinusoidal  components 
(caused  by  the  rigid  body).  Assume  that  the  indication  function  for  the  wide  win¬ 
dow  is  calculated  as  a  region  with  the  smallest  signal  samples  different  from  0.  We 
obtain: 
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However,  the  results  obtained  with  the  wide  window  are  not  a  good  indication  for 
a  region  of  rigid  body  part  since  they  include  the  influence  of  the  sinusoidal  FM 
signal. 
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Time-frequency  plane"  for  a  narrow  window 
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Average  for  the  largest  five  samples 
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Now  we  can  recall  three  functions  LI  obtained  for  three  various  windows. 
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To  obtain  the  detection  function  for  the  signal  region  caused  by  the  rotating  part 
it  is  enough  to  reconsider  the  time-frequency  representation  with  middle  window 
size,  i.e.,  to  consider  the  statistics  of  the  highest  and  median  samples  for  the  middle 
window. 


L, 

3V  uj  )  oicu 
0  1  2 

IMlt 

6 

3  1U1 

2 

lilt 

0 

1  CU  p  vjL 

3  |  5 

imp  1 

5 

to 

5 

7 

11 

7 

5 

5 

5 

3 

0 

0 

2  |  6 

2 

L; 

l 

L. 

2(03)  Median  of  samples 

0 

2 

6 

2  I  0  I  0  '  0 

0 

^0 

2 

6 

2  0 

0 

0 

0 

0 

0 

2 

6  1 

2j 

3(a))/L2((o) 

0/0 

1 

1  1  1  0/0 

inf  inf 

inf 

inf 

3.5 

;  i.83 

3.5  |  inf  inf 

*  inf  inf 

0/0 

0/0 

1 

II 

r’(w)=L3(co)/L2(co)>=(  l+p),  (5=  l  (0/0  treated  as  0) 

0 

0  0 

0 

0 

1  |  1  |  1 

l 

1 

0 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

The  pattern  caused  by  the  rotating  part  can  be  assumed  to  be  the  region  of  the  TF 
plane  distinct  from  0  with  I"{uj)  =  1 


DRDC  Ottawa  TM  2006-187 


33 


!’(©) 


l”(co) 


0 

0 

0 

0 

« 

t 

« 

1 

n 

I 

M 

1 

I  — 

1 

0 

l  ™ 

1 

H 

1 

s- 

L 

» 

i 

sm 

i 

0 

0 

0 

0 

J 

0 

blue  color  rigid  body  influenced  signal 

red  color  rotating  body  influence 


34 


DRDC  Ottawa  TM  2006-187 


The  following  Figures  A1-A13  illustrate  the  step  by  step  procedure  to  separate  the 
rigid  and  rotating  parts  of  the  body  using  the  method  described  in  Section  3.  A 
similar  signal,  as  in  Section  3,  is  used  in  this  example.  Figure  A1  shows  the  time- 
frequency  representation  of  the  signal  with  a  very  wide  window  and  Figure  A2 
shows  the  sorted  time-frequency  representation.  The  statistics  for  the  wide  window 
case  is  shown  in  Figure  A3.  Figure  A4  shows  the  time-frequency  representation 
of  the  signal  with  middle  window  and  Figure  A5  shows  the  sorted  time-frequency 
representation.  The  statistics  for  the  middle  window  case  is  shown  in  Figure  A6. 
Figure  A7  shows  the  time-frequency  representation  of  the  signal  with  narrow  win¬ 
dow  and  Figure  A8  shows  the  sorted  time-frequency  representation.  The  statistics 
for  the  narrow  window  case  is  shown  in  Figure  A9.  The  indicator  function  for  the 
stationary  patterns  (rigid  body  influence)  is  shown  in  Figure  10.  Figure  1 1  shows 
the  indicator  function  for  sinusoidal  FM  signal  patterns  (rotating  part  influence). 
The  reconstructed  rigid  body  and  rotating  parts  are  shown  in  Figure  12  and  Figure 
13,  respectively. 
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Figure  A. 1:  Time-frequency  representation  with  very  wide  window 
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Figure  A. 2:  "Sorted”  time-frequency  representation 
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Figure  A. 3:  Statistics  for  the  wide  window  case  (LI  smallest  samples ,  L2  median 
samples ,  and  L3  highest  samples) 
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Figure  A.4:  Time-frequency  representation  with  middle  window 
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Figure  A.6:  Statistics  for  the  middle  window  case  (LI  smallest  samples,  L2  median 
samples ,  and  L3  highest  samples) 
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Figure  A.7:  Time-frequency  representation  with  narrow  window 
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Figure  A.9:  Statistics  for  the  narrow  window  case  (LI  smallest  samples,  L2  median 
samples,  and  L3  highest  samples) 
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Figure  A.11:  Indicator  function  for  the  sinusoidal  FM  signal  pattern  (rotating  part 
influence) 


DRDC  Ottawa  TM  2006-187 


43 


500 

450 

400 

350 


Figure  A. 12:  Detected  stationary  pattern  (rigid  body  influence)  after  the  proposed 
procedure 
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Figure  A. 13:  Detected  sinusoidal  FM  signal  (rotating  part  influence)  after  the 
proposed  procedure 
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